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ABSTRACT 

The problem of eliminating the statistical fluctuations and extracting the event dy- 
namics from event-by-event analysis is discussed. New moments G p (for continuous dis- 
tribution), and G q ^p (for anomalous distribution) are proposed, which are experimen- 
tally measurable and can eliminate the Poissonian type statistical fluctuations to recover 
the dynamical moments C p and C qtP . In this way, the dynamical distribution of the 
event-averaged transverse momentum p t can be extracted, and the anomalous scaling of 
dynamical distribution, if exists, can be recovered, through event-by-event analysis of 
experimental data. 
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1 Introduction 



In the conventional investigation of high energy multiparticle production, all the events in 
a collision process are taken as a whole and the distributions, fluctuations and correlations 
inside the sample are studied without distinguishing the individual events. This kind of 
study is usually referred to as inclusive one. In the 80's - 90's of last century, motivated by 
the experimental observation [I] [2] and theoretical perspective |3] , people started to carry 
on the study event by event [I] [0] [Z] , with the aim of exploring the possible existence 
of new physics, such as quark-gluon plasma (qgp) [8] or non-linear dynamics 

One of the differences in these two kinds of study is that, the inclusive dynamics 
can readily be extracted from the experimental measurement through averaging over a 
large number of events, while the extracting of event dynamics is complicated due to the 
statistical fluctuations coming from the limited number of particles in a single event. 

For concreteness, let us consider the event-by-event fluctuation of the transverse mo- 
mentum distribution p(pt), a knowledge of which is important for the understanding of 
the basic collision dynamics [12]. Usually, it is convenient to "coarse-grain" p{p%), i.e. to 
divide the phase space region A of p t into M bins and integrate p(p t ) over p t in the mth 
bin 5 m 

p m = I p(pt)dp t , (m = 1,2,...,M). (1) 

The set 

p h p 2 ,...,p M (2) 

is the "coarse-grained" distribution 19]. When M — > oo (S m — > 0) it recovers the orig- 
inal distribution p(pt) [TT] . If p(pt) is continuous this process is convergent and a good 
approximation for p(p t ) could be obtained for not very large M. 

The realization of p(pt) in experiment is the distribution of the total number iV of 
particles in the p t region A, and the expression 

q m = NJN (m — 1,2, ... , M) (3) 

is used to evaluate p m , where N m is the number of particle falling into the mth bin. This 
is, however, exact only when iV — > oo. Just at this point appears the difference between 
the inclusive and event-by-event studies. 
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In an inclusive study, what is under consideration is the transverse momentum dis- 
tribution Pinci(pt) in the event sample. In this case, the number N in Eq.(3) is the total 
number of particles in the whole sample and could be made arbitrarily large through 
increasing the number of events in the experiment and thus the effect of statistical fluctu- 
ation can be gotten rid of through averaging over a sufficiently large event sample. On the 
contrary, in an event-by-event study the distribution p(p t ) in consideration is the trans- 
verse momentum distribution in a single event and the number N in Eq.(3) is the number 
n of particles in the event, which is limited by energy conservation, and the statistical 
fluctuation inevitably comes in. 

Various methods have been proposed to eliminate the influence of statistical fluctuation 
and evaluate the dynamical ones. Most of them are based on the comparison of the 
measured fluctuation with the expectation of statistically independent particle emission. 
For example, in Ref. [S] the results from mixed events are considered as the baseline for 
the random distribution and the difference in the fluctuation from a random distribution 
defined as 

d = l^data — ^baseline (4) 

is taken as a measure of the dynamical fluctuation. In Eq.(4) 



>/<Pt 2 > ~ (Pt) 2 v „ 
(Pt) (Pt) 

p t is the mean transverse momentum in a single event and (• ■ ■) denotes the average over 
event sample. 

Alternatively, in Ref. ^3] the statistical variance of event mean p t , under the assump- 
tion of independent particle production, is estimated as 



2 _ ^Ptincl /«\ 

°p- t stat - / \ J V°) 

and the difference between the variances of pt obtained from data and "stat" is taken as 
the dynamical variance 

2 

2 _ 2 _ 2_ __ 2_ _ ^Ptincl /~\ 

ptdynam ^ptdata pistat ptdata / \ ' \ ) 

A widely used measure for the non-statistical mean p t fluctuation is the $ Pt proposed 
in Ref. [12] 

$ pt = y/(Z*)/(n) - v^) , (8) 
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where z and Z are defined as z = p t — (pt) for each particle and Z = YJ}=i z % = n (pt — (pt)) 
for each event, respectively ^3]- The second term of the r.h.s. of Eq.(8) is the square root 
of the inclusive variance <Jp tincl = ((p t — (p t )) 2 ). Assuming that the multiplicity fluctuation 
is uncorrelated with the p t fluctuation, we get from Eq.(8) 



This equation is evidently similar to Eq.(7), both have the same structure as Eq.(4), being 
based on a subtraction procedure, i.e. to subtract the variance of p t or a quantity related 
to it, that will be expected from pure statistical fluctuation, from the same quantity 
obtained in experiment. These measures will, of course, vanish for a pure statistical 
system, and a non- vanishing value of them will indicate the existence of dynamical effect. 
Therefore, the measures based on the subtraction procedure, as those listed above, will 
at least qualitatively measures the effect of dynamical fluctuation. 

The aim of the present paper is to develop a systematic method, which is able to 
eliminate the statistical fluctuations directly from the experimental event-by-event anal- 
ysis and extract quantitatively the dynamical p t moment of any positive integer order, 
under the assumption that the statistical fluctuations are of the Poisson type, i.e. due to 
uncorrelated random particle emission. 

In Section II a theorem will be proved which is the basis of the elimination of Poissonian 
statistical fluctuations. In Section III the theorem is applied to the event-by-event analysis 
of transverse momentum distribution. In Section IV the event-by-event fluctuation of the 
non-linear fractal property will be discussed and the proposed method will be applied 
to this case to extract the dynamical fluctuation of fractal property. Section V is the 
conclusions. 



2 The elimination of Poissonian statistical fluctua- 



Divide the transverse momentum region A into M bins. Let p m be the event dynamical 
probability of pt in the mth bin, cf. Eq.(l), and n m the number of particle in the event 




(9) 



tion 
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with p t lying in the mth bin, then we have the theorem: 



/ Vf^\-/Vf W ™( W ™ l)---(n m p+1) i 



where f m is an arbitrary variable depending on m. 

Before going on to prove this theorem, let us notice that the symbol (■ ■ •) in the 
two sides of Eq.(lO), despite of both being the average over event sample, have different 
meanings. In the l.h.s. it is simply the average over dynamical probability distribution, 
while in the r.h.s. it includes also the average over Poisson distribution of particle number. 

I M \ r M 

( £ fmfra) = I £ /rn^P«d Pm , (11) 
\m=l / m=l 

/ V M f n m (n m -l)-{n m -p+l) \ 

r M 00 n (n —l)---(n -o+l) e -vm{n) 

= / E E fm m[ m \ m P } ■ -^(Pm(n)r-P(p m )d Pm , (12) 

J ra=ln m =0 \ n / n m- 

where P(p m ) is the dynamical probability distribution of p m in the event space, p m (n) = 
(n m ) is the average multiplicity in the mth bin. Using the normalization condition of 
Poisson distribution, it is easy to see that the r.h.s. of the above two equations are equal 
and Eq.(lO) follows. 

Note that a simplified version of the above theorem with f m = 1 in Eq.(lO) has been 
proved in Ref. P and has been widely used in the intermittency study [4 . Eq.(lO) is a 
more general theorem with a factor f m included, which is essential in the application of 
this theorem to the elimination of statistical fluctuation in event-by-event analysis. 

3 The elimination of statistical fluctuation in event- 
by-event transverse momentum fluctuation 

Let p t and pt exp be, respectively, the dynamical and experimentally measured values of 
the event-averaged p t , 

. m 

Pt= p t p(pt)dpt= £(pt)mPm, (13) 
JA m=l 



M M 

1 1 I, 



Ptcxp = E (Pt)m5m = E (Pt)m-T 1 , (14) 



m=l m=l 



where (pt) m is the p t value in the mth bin. The event-space moments of p t and pt exp are 

CM) = (Pt p ) = ( ( J2(pt) m pn)j y (15) 

C; XP (Pt) = ((plexpF) = ((l>t) m ^) ^ , (16) 

respectively. 

Let us first consider the elimination of statistical fluctuations in the second and third 
order moments. When p = 2, 3 we have 

= ( S (prfmPm) + ( E (Pt) m (Pt) m 'PmPm'\ , (17) 
\m=l / \m^m' I 



M 



C 3{Pt) = \\H (Pt)mP 



m 

\m=l / I 

M \ I M 



= ( 2 (Pt)mPm ) + 3 ( E (Pt)m(Pt) 2 m >PmP 2 m > ) 
\m=l I \m^m' I 

+ { (Pt)m(Pt)m'(Pt) m"PmPm'Pm" 

(18) 

\mj^m'^m" I 



Define 



G 2 (ft) = (|(ft) ^ m( "; )2 " ^ ) + (jtjP^Mrn^) , (19) 

G 3 (Pt) = (|>); ^ m(nm =^ ' 2) ) + 3 ( f ^ )m(Pt) ^ """"(iy " ^ ) 

+ ( e bt) m ( Pt )^(A) m » nm ;;: 3 nm " )- (20) 

It is easy to see, using Eq.(10), that G 2 (pt) = ^(pt), ^(pt) = ^(pt), provided the sta- 
tistical fluctuations are Poissonian. Thus the dynamical moments C2, C3 can be extracted 
from the experimental measurement by using G 2 , C3 and the statistical fluctuations have 
been eliminated. The elimination of statistical fluctuations in higher order moments can 
be proceeded in a similar manner. 
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Figure 1: Variance a^ t of p t distribution 



7 



Let us demonstrate the above results with a toy-model Monte Carlo simulation. In 
this model the distribution of pt is taken as 

P(pt) = ^Pte~ 2pt/a , (21) 
a 

with a Gaussian distributed parameter a (cx 2 (a) = 0.24). In total 1,000,000 events have 
been generated. The resulting variance a^ t = (pt 2 ) — (pt) 2 for C, C* exp and G are plotted 
in Fig.l as upward triangles, full circles and downward triangles, respectively. It can be 
seen from the figure that the width of experimentally measured p t distribution is 
wider than that of the dynamical distribution, especially when the average multiplicity is 
low, while the width calculated from G coincides with the dynamical width, er|~ = er^ , 
where a^f = G 2 (p t ) - (Gi(p t )) 2 and Gi(p t ) = (E TO (Pt)mn m )/(n). Therefore, using G p 
the Poissonian statistical fluctuations are eliminated thoroughly and the event dynamics 
is successfully extracted. 

4 The elimination of statistical fluctuation in event- 
by-event analysis of fractal property 

In the above discussion the dynamical event-distributions of phase space variables - 
rapidity y, transverse momentum p t , azimuthal angle ip - - are implicitly assumed to be 
continuous functions, i.e. fulfil the condition 

5x^0 =^ Sp(x) -> 0, x = y,p t ,(p, (22) 

which means that when the bin size decreases the probabilities in neighboring bins tend 
to be equal to each other. 

However, there is evidences showing that this may not be true in some cases. The first 
experimental evidence is from a jacee event in 1983 jT| , in which the total multiplicity 
is about one thousand and the multiplicity fluctuations in a small rapidity bin are still 
2 x the average. Similar phenomena have also been observed afterwards in accelerator 
experiments with local fluctuations up to 60 x the average j2]. Obviously, this kind of 
fluctuations is out of the usual statistical ones, indicating the existence of self-similar 
fractal property in phase space distribution. 
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A characteristic phenomenon of self-similar fractal is the anomalous scaling of the 
normalized probability moments (npm): 

C q {M) = M _1 / ^2 {M Pm ) q ) oc M+*. (23) 

\m=l / 

Such an anomalous scaling property can be demonstrated using a toy model for self-similar 
fractal - - random cascading a model [0] JH] . This model describes each multiparticle 
event as a series of steps, in which the initial phase space region A is repeatedly partitioned 
into A = 2 parts. After v steps we get M = 2 U sub-cells of size 5 = A/M. At each step v 
the probabilty in each of the two parts is obtained by multiplication of the probability in 
the step v — 1 by a particular value of the random variable wj^, where j v is the position 
of a sub-cell at the vtW step (1 < j v < 2 U ). The elementary fluctuation probability u can 
be chosen in various ways. The simplest way is to choose it as jT3] 

"$.1 = ^(1 + or), Jg = \{l-<*rl (24) 

where r is a uniformly-distributed random number in the interval [—1, 1], j is an integer 
(1 < j < 1 v ~ l \ a is a characteristic parameter of the model taking value in the interval 
[0,1]. 

In this model, after v steps of partition, the npm defined in Eq.(23) becomes 

C q (M) = M q (u q (l)---u q (v)). (25) 

Since (u) = 1/2, we have the anomalous scaling of NPM, cf. Eq.(23), 

C q {M) = M q {uj q Y = M q e uln{ " q) = M^, (26) 

with M = 2 U = e uln2 ,(p q = q + \n(u; q ) / In 2. 

It is easy to see using Eq.(lO) that the npm C q can be extracted from experimental 
data through the normalized factorial moments (nfm) 

F q (M) = M q ~ x f ^{n m -l)---{n m - q + l)) 

The anomalous scaling of nfm, usually referred to as intermittency |H]j3], has been ob- 
served in hadron-hadron jlH] and e + e~ ^H] experiments, providing a first signal for 
the non-linear fractal property of strong interaction dynamics. 



The npm defined in Eq. (23) can be viewed as the sample average of "event-probability 
moment" (epm) 

m 

Cj 1 c \M) = M- 1 J2(M Pm y. (28) 

m=l 

It is natural to consider the distribution of epm itself instead of only its average fl9j. 
This distribution can be characterized by the pth order event-space moment C qtP of the 
gth order event-moment C^> 

C q>P (M) = ((^ C) (M)) P ) = / (M- 1 £ (Mp m )<^ \ (29) 
and the corresponding normalized moment 

c£r ] = c^/tc^y. (so) 

In order to measure C q ^ p in real experiments, the pth order event-space moment of the 
gth order event-factorial-moment 

n,P = ( (m- t (3D 

has been proposed [19J. The scaling property of F q ^ p with the increase of M is referred to 
as erraticity [19J and has been observed in experiments [20J[21J. However, when p ^ 1, 
F q>p contains statistical fluctuations j22] and is unequal to C q>p , despite of the factorial 
moments apparently used. The observed erraticity phenomena have been shown to be 
dominated by statistical fluctuations and the dynamical effect is hidden [22 J [23 . 

In order to be able to observe the dynamical "erraticity", the statistical fluctuation 
should be eliminated first. The elimination of (Poissonian) statistical fluctuation for the 
case p — 1 is straightforward using Eq.(10). The result is just Eq.(31) with p — 1, or 
Eq.(27). Let us consider the elimination of statistical fluctuations for the cases q = 2, 
p — 2, 3. We have 

//A/\ 2 \ I M M \ 

C 2 ,2 = ( [MY.pI ))=M 2 (Y,Ptn+ Ett > (32) 

\ \ m=l / / \m=l ra^-m! I 

I I M \ 3 \ I M M M \ 

^2,3 = ( M Y. Vl ) = M 3 ( £ pi + 3 £ P 2 mPtn> + £ pIpIpI" ) • (33) 

\ \ m=l / / \m=l m^m! rn^m'^m" I 



Define 



j — i n/r 2 I \ " n m ' ' ' ( n m 3) \ 2 / x - n m (n m l)n m '(n m 1) \ 

2 = M \S w / + M \JL w / ■ im 

in 



3 = M L?, w / +3M \J& / 

, #3 / ^ n m (n m - l)n m /(n m / - l)n m »(n m » -1) \ 

\ A „ w /' 

Then, utilizing Eq.(lO), it is ready to show that G 2j2 = C 2i 2, G 2i 3 = C 2i3 , provided the 
statistical fluctuations are Poissonian. This means that using G 2j2 , G 2j 3 instead of F 2 2 , F 2 3 
to measure the dynamical moments C 2j2 , C 2j 3, the (Poissonian) statistical fluctuations are 
eliminated and the real dynamical "erraticity" is observed. This method can be extended 
to any positive integer orders q and p. 

To illustrate the above result by Monte Carlo simulation, the random cascading a 
model described above is used with a Gaussian distributed parameter a (mean=0.5, 
width=0.15). Particles are put into each event according to Poisson distribution with 
(n) = 6. The results from 2,500,000 generated events are plotted in Fig. 2. The coinci- 
dence of Gg°p m and C°° rm is remarkable. It is to be contrasted with the strongly upward 
bending curves of lnFq° rm vs. InM. This shows that the strongly upward-bending be- 
havior of lni^° rm vs. InM for positive p which is typical in the conventional erraticity 
analysis 19 20jj2I] is due to statistical fluctuations. The reason is: when the partition 
number M increases, the average multiplicity per bin becomes smaller and smaller and 
the statistical fluctuations get stronger and stronger. After eliminating the statistical 
fluctuations, the lnG^ p rm vs. InM curve coincides with lnC^° rm vs. InM and recovers the 
anomalous scaling property of the event dynamics (dynamical erraticity) in the model. 



V Conclusions 

In this paper the problem of eliminating the random noise in event-by-event analysis 
is considered. New moments G p (for continuous distribution), and G qyP (for anomalous 
distribution) are proposed, which are experimentally measurable and can be used to 
eliminate the Poissonian type statistical fluctuations and recover the dynamical moments 
Cp and Gqp. 

For a comparison of different methods we notice that most of the measures proposed 
in the market for the dynamical fluctuation of transverse momentum are based on a 
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Figure 2: Scaling property of event-space moments 

subtraction procedure, i.e. to subtract the variance of p t or a quantity related to it, that 
will be expected from pure statistical fluctuation, from the same quantity obtained in 
experiment. In the present paper we take another approach, i.e. not to evaluate the effect 
of pure statistical fluctuation and subtract it, but to eliminate the statistical fluctuations 
directly and quantitatively from the experimental data and recover the dynamical a^ t 
under the unique assumption that the statistical fluctuation is Poissonian, i.e. due to 
uncorrelated random particle emission. 

The difference between these two approaches is two-fold. 

Firstly, in the subtraction method how to get the "pure statistical" variance from ex- 
perimental data is a big problem. The mixing-event method [H] can in principle get the 
statistical variance, but the accuracy depends on the mixing-procedure. The method pro- 
posed in Ref. [Hj relies on the equality cXp- stat = &p tincl /(n), cf. Eq.(6), but this equality 
holds only for a pure statistical system without any dynamical fluctuation. In experimen- 
tal data sample there exist simultaneously dynamical and statistical fluctuations, and the 
statistical variance cr?- d s ^ included in the data sample is unequal to the inclusive variance 



I 1 ? 



of the data sample over average multiplicity: 

2data / 2data //_\ loa\ 

The same holds also for the $ Pt method [P2j . 

Secondly, the variance cr^. contains the p t moment only up to second order, which only 
gives the width of the distribution of p t in the event space. Using the method proposed in 
the present paper the pt moment of any integer order can be extracted and the dynamical 
distribution of p t in the event space is portrayed in much more detail. 

In the anomalous case, the difficulty coming from the dominance of statistical fluctu- 
ation in the theoretical and experimental studies of erraticity up to now can be overcome 
and the dynamical erraticity, if exists, can be extracted using the method proposed in the 
present paper. 

Applying the proposed method to real experimental data is highly recommended. 
Acknowledgement The authors thank Wu Yuanfang and Liu Feng for helpful dis- 
cussions and comments. 
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